Collisions of bosonic ultracold polar molecules in microwave traps. 

Alexander V. Avdeenkov 1,2 
1. National Institute for Theoretical Physics (NITheP), 
Stellenbosch Institute of Advanced Study, Stellenbosch 7600, South Africa 
2. Institute of Nuclear Physics, Moscow State University, 
Vorob'evy Gory, Moscow, 119992 Russia 

Abstract 

The collisions between linear polar molecules, trapped in a microwave field with circular polar- 
ization, are theoretically analyzed. The microwave trap suggested by DeMille l| seems to be rather 
advantageous in comparison with other traps. Here we have demonstrated that the microwave trap 
can provide a successful evaporative cooling for polar molecules in a rather broad range of frequen- 
cies of the AC-field. We suggested that not only ground state polar molecules but also molecules 
in some other states can be safely trapped. But the state in which molecules can be safely loaded 
and trapped depends on the frequency of the AC-field. 
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I. INTRODUCTION 



In the last few years a considerable progress has been made in the development and inves- 
tigation of various types of electromagnetic traps for different but especially polar molecules. 
The building of such traps should have made it possible to perform many interesting experi- 
mental studies and applications involving the manipulation of cooled and trapped molecules 
with kinetic temperatures near zero Kelvin Schemes for creating 

samples of cold molecules are being developed now in several laboratories around the world 
where simple molecules at kinetic temperatures in the mK to fiK range are being created 
and studied. One of the main purposes in the perfection of traps is to keep samples of 
molecules dense and cold enough and the microwave trap seems to be very advantageous in 
restraining polar molecules 1|. 

In this article we are studying the collisional dynamics of linear polar molecules in the X S 
state trapped in a microwave field with a circular polarization at cold temperatures. The 
choice of linear polarization does not seem very appealing as it was shown in [l[ as AC- Stark 
levels have a lot of avoiding crossings which enhances the collisional loss. Polar molecules 
have strong permanent electric dipole moments and therefore the collisional dynamics is 
mostly ruled by the dipole-dipole interaction: 

fiife - 3(e R /li)(e R /I 2 ) 



i? 3 



(1) 



Here, /2j is the electric dipole moment of the molecule i, R is the intermolecular separation 
and en = R/R. Thus in our model collisions are controlled by two ratios vjB and x = 
liE/hB (y is the microwave frequency, B is the molecular rotational constant, and £ is 
the electric field strength). We are using the coupled-channels method to calculate elastic 
and inelastic cross sections in dependence on the frequency and the strength of an AC 
electric field. It has to be noted that collisional dynamics of cold polar molecules in the 
DC- electric field was rather intensively theoretically explored which revealed quite rich 
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One of the main goals is to understand if polar molecules in their absolute ground state can 
have a successful evaporative cooling. At first sight it seems rather an inappropriate question 
as molecules are in their absolute ground state. Naturally, if the energy of the first excited 
rotational state is smaller than the microwave frequency(z/ < B) it seems that the system is 
safe against two-body inelastic collisions and collision losses are excluded. And this is true 



for bare particles. But it might be not completely true for particles dressed with a microwave 
field , especially for parameters which provide a large and deep trap(z/ ~ B,fiS <^ B). As 
we have seen the reason is that the AC- field and the dipole-dipole interaction are mixing 
states with different angular moments and "dressing" . We have found that even for the most 
favourable case v < B, there are inelastic processes. But these processes are causing mostly 
the "undressing" or "dressing" of the colliding particles without changing their internal 
structure. Here we have to be careful about the term "mostly" as dressed states are the 
superposition of lots of pure states. But what we can be sure of is that these " dressing" - 
"undressing" processes do not cause the loss of molecules from the trap. 

Moreover this is supposedly be true not only for molecules in their absolute ground state 
| J = 0, M = >, but for example in the \J — 1, M — —1 > state at frequencies around 
vjB = 3, for molecules \ J = 2, M = —2 > at frequencies around vjB = 5 and so on(J- 
rotational quantum number). Actually how far this "around" goes depends on the strength 
of the AC- field. 

We have found that during the collisions the parity of the projection of the total mo- 
mentum of two colliding dressed molecules is conserved. Our calculations demonstrate that 
the cross sections for such positive and negative projection parities are different by several 
orders of magnitude. 



II. POLAR X S- TYPE MOLECULES IN A MICROWAVE FIELD 

Here we analyze the AC-Stark shift, which determines the trap depth. The energy levels of 
1 £- type molecules can be described by rotation J, total spin F (i.e., including nuclear spin), 
and vibration v quantum numbers. For simplicity we will neglect hyperfine splitting and 
consider molecules only in the v = vibrational ground state. So we treat polar molecules as 
rigid rotors with a permanent dipole moment. The dressed- state formalism is particularly 



20] where 



convenient for describing the photon-atom interaction in the strong field limit 
the dressed states are eigenstates of the Hamiltonian of the total system: particle plus 
photon. The basis states of the rigid rotor plus field Hamiltonian, H = H rot + Hf ie i d 



are 



J, M, N + nj, where M is the projection of J on space-fixed axis, which is conveniently 
chosen parallel to the wave- vector of the microwave field. , N + n is the dressed state photon 
number and n << N is the deviation of the photon number. The AC- Stark splitting is 
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caused by the molecule- field interaction H S tark = —p£- Considering the circularly polarized 
microwave field (cr~), the non-zero Hamiltonian matrix elements (normalized by HB) are 
given by 

(J, M, n \H rot + H field \ J, M, n) = J ( J + 1) + (2) 



(J + l,M + l,n+ 1 LHsterfc J,M,n) = , = 



(J+ 1,M - l,n- 1 M, n) = £2^™^1, (3) 

Here, it is assumed that n is much smaller than the mean photon number N. In the AC- 
electric field the J, M, N + n are not good quantum numbers and the dressed state formalism 
should be applied. 

A. Dressed states 

The AC- electric field mixes states with different J, M and n and hence neither of them 
is a good quantum number. Thus we can only mark our states by their origin at zero 
field(Fig{T]) where they can be assigned by J and M quantum numbers. Figure Q] shows the 
energies of dressed states versus applied electric field frequency for the J = and J = 1 
states at x = 0.7 within v/B — [0, 4] frequency range. 

First we would like to consider the ground and so the lowest- field seeking state \J = 
0, M = >. Though the term lowest-field seeking state does not seem appropriate in dressed 
state formalism. In practice we transform the molecular state to a field-dressed basis for 
performing scattering calculations and the state with a given mean photon number N is 
described as: 



(JMn)) = (JMn);£) = £ p (JMn) \JMn) , (4) 

JMn 

where P(jMn) stands for the eigenfunctions of the H rot + Hf idd + H stark Hamiltonian de- 
termined numerically at each value of the field. We will continue to refer to the molecular 
states by the quantum numbers J, M, and n, with the understanding that they are only 
approximately good in a field, and that Eq. [4] is the appropriate molecular state. Note that 
AC- Stark splitting should be irrelevant to the value of |n| as n << N so we choose small 
values of n. Further we are mostly interested in cases in which the strength of the AC- field 
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FIG. 1: Energies of dressed states versus applied microwave electric field frequency at x = 0.7. 
See text for details. 

is not large(/i£ < B) which is found to be more practical for a microwave trap [li]. In this 
case the above mentioned mixing is not that strong and allows us to analyze different cases. 



B. 



< 2B 



case 



(000)) state which is a pure |000) if we neglect the molecule- field 



Let us consider 

interaction H stark and describe molecules in their absolute ground state. In [l| it was 
pointed out that the most relevant trap parameters are at x ^ 1 and the ratio u/B ~ 1. In 
this case the state of interest is approximately described as: 



(000))^p ( ooo)|000)+p ( 



l-i-i) 



1) 



(5) 



It has to be noted that at x ^ 1 , P(ooo) ~ P(i-i-i) ~ l/\/2. At larger x admixtures of other 
states, like |2 — 2 — 2) and so on are not small. Here we would like to point out that if trapped 
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molecules are in this state (000) ) it is equal to saying that trapped molecules are in (00|n|) 
states as the AC Stark shift is the same for all of them. Moreover the H rot + Hf ie id + H stark 
Hamiltonian does not mix, for example, states like (000)) and (00 ± 1)). But the dipole- 
dipole interaction can mix them as the 



(000)) (000))^ (00 ±1)) (00 ±1) 

scattering processes are possible. In order to reveal the possibility of such processes one has 
to see that at x ^ 1 and the ratio vjB^X these states are mostly described as: 



(oo - r 



P(oo-i) |00 - 1) + P(i_i_2) 11-1-2) 
(001)) wp (00 i)|001)+p ( i_io) |1-10) 



Thus the 



(000) fi (00 - 1)) ~ P(i-i-i)P ( oo-i) (1 - 1 - 1| fi 1 00 - 1) 
(000) n (001)) ~ p ( ooo)P(i-io) (000|// |1-10) 



(6) 
(7) 

(8) 



matrix elements provide for the above-mentioned process and their values depend on p^jMn) 
eigenstates. 

At v % B the only two inelastic processes may have considerable cross sections: 

(9) 



(000)) (000)) => (00 - 1)) (001) 



(000)) (000) W (00 - 1)) (00 - 1) 



Note that the threshold energies of (000)) (000)) and (00 - 1)) (001)) channels are de- 
generate. 

At v > B the 



1 — 1 — 2)) state becomes energetically accessible (see Fig. CD) so that the 
(000)) (000) W (00~1)) (l-T-2)) (10) 



inelastic process is possible. 



C. 2B<v<AB case 



In this region of frequencies our state of interest (000)) is approximately described as: 



(000)) w p ( ooo) |000) |1-1-1) +P( 2 -2-2) |2-2-2) 



(11) 



If the strength of the field is not large(see Fig. [T]) then in the middle of this region this state 
is an almost pure |1 — 1 — 1) state. As vjB gets close to five, the state is dominated by the 
1 2, —2, —2) state. So the larger the frequency of the AC- electric field, the larger particular 
angular moment J will be dominated. But this is only true if the strength of a field is not 
large(x < 1). Otherwise, at a large x, the dressed state will be a superposition of a lot of 
states. 

Compared to the v < IB case there is one extra open channel for inelastic collisions: 

(000)) I (000)} |(1 -T- 2)) |(1 -T- 2)) . (12) 

Moreover our state in this region of frequencies becomes a weak-field seeking one while the 
state (1 — 1 — 2)) becomes a strong-field seeking state. 



III. RESULTS AND DISCUSSIONS 

The main difference from the DC electric field case is that in general inelastic processes 
are always allowed for any non-zero field and for any state. Of course an inelastic cross 
section will be very small at small fields but not exactly zero. In the following we focus only 
on bosonic species. The case of fermionic species will be considered elsewhere. 

We have solved the coupled- channel equations using the logarithmic derivative prop- 



agator method [21J to calculate total state-to-state cross sections. The scattering calcula- 
tions quickly become computationally expensive as more rotational states, partial waves 
and dressed-state photon numbers are included. We have found that it is enough to take 
J = 0, 1; I = 0, 2 for bosons and I = 1,3 for fermions and An = 0,±1,2 for changing 
photon numbers to get reasonably converged cross sections. The only two-body interaction 
taken into account here is the dipole-dipole interaction which is sufficient to describe the 
qualitative behavior of cross sections. 

We have not chosen any particular polar molecule for calculation but a rigid rotor with 
the rotational constant B = 6GHz and dipole moment \x = 3D. For this molecule, we are 
mostly interested in the lowest energy strong-field-seeking state of the ground vibrational 
state. The cross sections were calculated for dressed molecules, which means we show them 
for molecules in (JMn); \ basis. The trap depth for this particular molecule is YlmK at 
x = 0.5 and 50mA" at x — 1 if vjB — 1/3. It is 60mA at x — 1 and vjB — 1. For the high 
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FIG. 2: Elastic(solid lines) and inelastic (not solid lines) cross sections versus the collisional energy 
for x « 1, v/B ~ 1. The black and red curves are for the positive and negative projection parities 
respectively. The dashed curves are for (000)) (000)) =^ (00 — 1)) (001) ) , the dotted curve is 



for 



(000)) (000) 



(00 - 1) ) (00 - 1) ) and the dot-dot-dashed curve is for (000) ) (000) 



(00 - l)j (1 - 1 - 2)j transitions. 

frequency case v/B = 2.5 it is only 3.5mK at x = 1. 

We have found that if one takes into account only the dipole-dipole interaction in the 
presence of an AC-field the only conserved value for two molecules in (JiMirii); (t/2^2^2)) 
states is (_i) M i+ M 2+m;+m+n2 -which we call a projection parity. Thus the system of the 
coupled- channel equations is split into two blocks with a positive and negative projection 
parity. The nature of this symmetry will be described elsewhere. 

Fig. [2] shows cross sections for bosonic and fermionic species for the (000); ) state at 
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x = 1 and u/B = 1 parameters which are proposed as more relevant for the future ex- 

n 

periment The black and red curves are for positive and negative projection parities 
respectively. The largest inelastic cross section is for (000)) (000)) =>- (00 — 1)) (001)) 
transitions and moreover it is not small and the threshold behavior is similar to that of an 
elastic process. The latter is clear as these states are degenerate. Such inelastic transition 
can be described as the "dressing" or "undressing" of molecules by one photon each. Will 
this "undressing" cause the loss of molecules from the trap? We believe that it will not as 
in the limit of a strong field the Stark shift as well as the trap depth will be the same for 
dressed states with N and N ± 1 photons. The only inelastic process which causes the trap 
loss is | (000)) | (000)) |(00~ 1)) |(1 -T- 2)). But the cross section for it is relatively 
small in the region from about 0.1/iK to IOO^lK and even maybe at lower temperatures. 
One can see that there exists some resonance structure. We do not take any notice of these 
resonances as they are defined by an unknown short-range part of the interaction which in 
our case is simply modeled by the inner wall of the interaction. Thus after all our "abso- 
lute ground state" molecules undergo inelastic transitions. Though at x ~ 1 and u/B pa 1 
trapped molecules will feel pretty safe. 

It is noteworthy to note that the elastic cross sections are much larger for positive parities 
at ultracold energies(around 1/J-K) especially while they are more or less of the same order 
at cold energies (around lmK). We have found that the threshold behavior of the inelastic 
cross section for bosonic molecules with a negative projection parity does not follow the 
Wigner law at the shown collisional energies at least. The only reason we can think of is 
that for a given projection parity not all projections mi of the partial wave I are involved. 
This phenomenon will be studied elsewhere. 

Fig. [3]shows cross sections for x = 0.5 and different u/B ratios . Though the trap depth is 
not that deep for this x it is still several mK. We would like to pay attention to the relatively 
small value of the field strength as in this case dressed states are mostly superpositions only 
of only two pure states or even sometimes one pure state (see discussions in Section II). In 
this case it is easy to analyze hidden details. At u/B pa 0.3 ratio the field is rather weak so 
dressed states are almost pure states( for example, it means that, say, (00 — 1)) ~ | (00 — 1)), 
see Eq.(5-8)) and all inelastic cross sections are small. At u/B pa 1 ratio the dressed states 
are mostly superpositions of two states(Eq.(5-8)) which considerably enhance "dressing"- 
"undressing" inelastic processes involving considerable transitions to |J = 1,M — — 1 > 
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FIG. 3: Elastic (solid lines) and inelastic (not solid lines) cross sections versus the collisional 
energy at x 0.5 for a)u/B = 0.3, b) 1 and c) 2.5. The dashed curves are for (000) ) (000) 



(00 - 1)) (001)) , the dotted curves are for (000)) (000)) =► (00 - 1)) (00 - 1)), the dot-dot 
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dashed curves are for 



(00 — 1)) (1 — 1 — 2)), and the dot-dashed curve is for 
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(1-1-2)^ (1 - 1 - 2)j transitions. 
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FIG. 4: Elastic(solid line) and total inelastic (dotted line) cross sections versus the microwave 
electric field strength at the collision energy lfiK for v/B = 0.3. The value x « 1 is attained at 
£ = AW /cm. 

states by means of the dipole-dipole interaction(Eq.8). Nevertheless such transitions are 
between equally trapped states (00|n|)^. At v/B = 2.5 ratio the dressed state of interest 
is by the largest part the |(1 — 1 — 1)) state and it is already a weak-field seeking state. 
Moreover one can see that inelastic cross sections which may cause the loss from a trap 
are really small. We may suppose that it is possible to prepare polar molecules in their 
| J = 1, M — — 1 > states and then load them into a microwave field at a frequency v around 
3B. It is probable that the same trick can be done with, say, J = 2, M = — 2 molecules by 
loading them into a microwave field at a frequency v around 5B. 

We have found out that the cross sections have rather a weak dependence on the frequency 
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FIG. 5: Elastic(solid line) and inelastic (not solid lines) cross sections versus the microwave electric 
field strength at the collision energy \\xK for u/B = 2.5. The inelastic processes are marked as on 
Fig. El 

of an AC- field (Fig H|5|) for large values of the strength of a field(x ~ 1 at E m AkV/cm). 
Both elastic and inelastic cross sections are more or less of the same order. Of course 
there is rather a rich resonance structure of cross sections which reflects the fact of the 
existence of lots of bound states in closed channels. We do not pay attention to this fact 
as the short-range potential is unknown and simply parameterized by the inner wall. But 
we would like to note that the positions of resonances will depend not on the strength of a 



field(like it is for DC- field [221 ]) but also on a frequency. Thus simply changing the frequency 
of a field one can catch a resonance. This issue will be under our investigation elsewhere. 
FigJH shows that at large values of field the inelastic "dressing" -"undressing" cross sections 
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become comparable with the elastic one but the low frequency (V/.B = 0.3) does not allow 
other trap losses processes. Thus the evaporative cooling may be successful for rather large 
strengths of AC-field as well. But it has to be checked as for a larger field we have to take 
a larger Hilbert space of basis states | JMn >. It means that some other inelastic channels 
with larger J,M and n may appear at a larger field and "destroy" such a pleasant picture. 
But it will be very costly numerically as we will have to handle up to several thousands of 
coupled-channels equations. Figj5] shows the same pattern of behavior for a higher frequency. 
The only difference is that there already are some inelastic channels at large fields which 
definitely cause trap losses. 

To recapitulate, we have computed scattering cross sections for cold collisions of the 
bosonic polar 1 £- type molecules as functions of both collision energy and AC electric field 
parameters. We found that even ground state molecules have inelastic transitions. Though 
most of them are harmless and allow keeping molecules in a trap. We suggested that 
polar molecules at \J, M = — J > states may have a successful evaporative cooling in the 
appropriate region of frequencies of a microwave field and will be restrained and shielded by 
this field against collision losses. 

We acknowledge useful interactions with Masatoshi Kajita and Masato Nakamura. 
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